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Abstract 

We study the passage of UHE neutrinos through the Earth in order to perform an absorption tomography of its inner 
structure. The aim of this work is to study the extraction methods of the Earth's density, in this conditions, we do not need 
to implement a realistic Monte Carlo simulation, as we are only interested in comparing the goodness of a standard method 
[1] with the one we propose. The Earth's density is reconstructed using the 2-d Radon transform and we compare the density 
obtained considering neutral current regeneration through the complete transport equation, with the one obtained making use 
. of the effective cross section approximation (standard method). We see that the effective cross section leads in general to 
' inaccurate results, especially for flat initial neutrino fluxes, while the full transport equation method works regardless of the 
^ , initial flux. Finally, an error propagation analysis made for different uncertainties in the surviving neutrino flux shows that the 
^ recovered density presents a percentage uncertainty less than two times the uncertainty in the flux. 
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Neutrino absorption tomography has been considered as an alternative way to obtain information about the interior 
of the Earth that would be independent of any geophysical model (see [1] and references therein). In view of this, it is 
useful to understand the role of neutral current regeneration of neutrinos in their passage though the Earth. To that 
\ end, in the present work we perform the absorption tomography following two different approaches to take account 
of the neutral current regeneration effect in the density extraction method. The first approach, which we call the 
O i' standard method or the effective cross section approach [2] , is largely used in the literature consists on implementing 
D . an approximate solution to the transport equation by defining an effective cross section. The other method makes 
use of the complete integro-differential transport equation itself and not of its solution, and this is the approach we 
introduce in this work. 



In the next section we start by describing both the standard and the transport equation approaches, and apply 
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H , them for some initial neutrino fluxes. We use theorized neutrino fluxes in the range lO^'GeV to lO^GeV coming from 
■ - - 1 cosmic sources such as Active Galactic Nuclei (AGN) and Gamma Ray Bursts (GRB) as well as neutrinos generated 
by cosmic ray interactions in the atmosphere (ATM). All these fluxes are assumed to be isotropic, and the Earth 
density is considered spherically symmetric so that the angular data of the surviving flux can be used to recover the 
density be means of the 2-d Radon transform which is a standard technique, but as we have not found it explained in 
any article, we decided to include it in the appendix. In order to have a surviving neutrino flux after traversing the 
Earth, we consider the numerical solution to the complete transport equation for neutrino propagation. This topic 
is presented in section 2 and the reconstruction procedures are discussed in section 3, where we compare the results 
of our approach with the ones of the (Te^j- approach. For the sake of comparing the two density extraction methods, 
it suffices to consider the numerical solution for the surviving flux, since a realistic Monte Carlo study would be 
consistent with the numerical results for the surviving flux. Finally, in section 4 the error propagation is performed 
when recovering the Earth's density following the transport equation approach, that is making no approximation 
related with the NC regeneration. 
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II. NEUTRINO PROPAGATION THROUGH THE EARTH. TRANSPORT EQUATION AND 

EFFECTIVE CROSS SECTION 



As neutrinos travel through the Earth they may suffer charged current (CC) and neutral current (NC) interactions 
with the miclcons in their path. Neutrino oscillation within the Earth can be neglected for energies higher than 
lOOOGeV [3], and since the CC and NC cross sections increase with the neutrino energy, we choose our energy range 
from lO^GeV to lO^GeV where the Earth is neither totally transparent nor opaque to neutrinos. 

The change in a neutrino flux (j)u{E) = as it traverses the Earth can be divided into two effects: absorption 
and regeneration [4] . Absorption is a decrease in the neutrino flux due to CC or NC interactions When neutrinos pass 
through an amount of matter dr = n{r)dz in a distance dz, where n(r) is the Earth's number density, the change in the 
flux 4>„{E, t) due only to absorption is proportional to (t>i^{E, r) by the total cross section atot{E) = acc{E) + (Jnc{E) 
which represents a probability of CC or NC interaction: 

= -atot[E)(j)^{E,T) (1) 

Here t{z) is the amount of material found until a depth z that is, 

t{z) = r dz'n{z') (2) 



Jo 

where the number density is the Avogadro's constant times the density, n(z') = Nap{z'). 

To write the complete transport equation for neutrinos of our energy range, we have to add to (1) the effect of 
regeneration, which accounts for the possibility that neutrinos of energies E' > E may end up with energy Ei, due 
to NC interactions with the micleons. Unlike absorption, regeneration represents an increase in the flux at an energy 
E caused by NC interactions at energies E' > E. Now, this change due to regeneration is proportional to each 
flux 4>iy{E' , t) by the probability of such interactions which is given by ^{E', E) so that we can write the transport 
equation as follows 



'{E.t) 



dr 



= -atot{E)ME,r) + dE' ' ME',r). 



(3) 



to be solved with the initial condition (j)v{E,Qi) = $,^(£'), where ^v{E) is the initial flux. 



The other approach consists on defining an effective cross section (jeff{E,^v) in an attempt to incorporate both 
absorption and regeneration all together. Dividing the transport equation (3) by (^^(S, r) we obtain 

d\n4>,{E,T) rp,^ r^p' iA\ 
- -<Jtot{E) + / dh — — -T-TTT— r- (4) 



dr ' Je^ dE ME,t) 

In this last expression, the quotient of the fluxes '^-j appears integrated on E' , and these fluxes must be evaluated 
at a depth in the Earth corresponding to r nucleons traversed. The effective cross section approach consists in replacing 
the above quotient by the initial fltix ratio [2] 



(5) 



and now we identify the effective cross section as 

t^^^ , r danc{E',E) <P,{E') 
a^ffiE, = -a,M + J^^ dE (6) 

The effective cross section results then from the approximation (5) to use the initial fluxes (which are supposed as 
known) rather than the fluxes (pv{E, r) at different depths through the Earth which can be obtained using the original 
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transport equation (3). When using the effective cross section, the resulting flux for an amount of material r traversed 
is 



ME,r) = ME)e-''^ff^''^- (7) 

We can compare this prediction for the surviving flux with the one obtained by solving the transport equation 
numerically. To do this, wc use the different initial diffuse fluxes for: ATM neutrinos derived by Volkova [5], neutrinos 
coming from AGN according to Stecker-Salomon [6] and to Protheroe [7], and neutrinos from GRB derived by Waxman 
[8] (see figure 1). 

Besides the initial flux, we also need to assume a density for the Earth to solve the transport equation or to use the 
(Te//-approach, so wc take this density to be given by the Preliminary Reference Earth Model (PREM) [9], which is 
plotted in figure 2. Neutrinos will then follow a path through the Earth towards the detector as illustrated in figure 
3. 

Solving the transport equation by Euler's method and by the iterative method described in [10] give the same 
results which are also consistent with the ones in [3]. For illustration we show in figures 4 and 5 the normalized final 
flux or shadowing factor 

as well as the shadowing factor Sgf / according to the (Te/ / approach 

SeffiE,e) = e-^'ff^^'^''^^ (9) 

for nadir angles 6 = {80°, 40*^, 0°} from the downward normal to the surface where the detector is (figure 3). 

We note that the two approaches agree when dealing with fast decreasing initial spectra such as ATM, while for 
flatter initial fluxes, the aef / approach trends to overestimate the surviving flux, (and hence the shadowing factor) 
as the nadir angle decreases. 
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III. TOMOGRAPHY OF THE EARTH 



In the appendix we see that the Radon transform of the number density is the amount of material traversed r 
appearing in the transport equation, and the expressions for the Radon transform and its inverse transform are 
obtained exploding the assumed spherical symmetry of the density. For reconstructing the number density we will 
then use (see the appendix) 

1 de' dT{e') 

n(r) = / =^—^ (10) 

Jarcsin ^ Vi?2 sin^ 0' - dO' 



This expression relates the number density of the Earth at a distance r from its center, with the different amounts of 
matter traversed at different nadir angles through the Earth. Then, measuring t(6') at all the possible nadir angles 

between y tt/2 wc may obtain n(r). 

To use equation (10), we must be able to infer ^ from the data of the surviving neutrino flux. We may then do 
this inference under the transport equation approach or under the c7e//-approach and reconstruct a density in each 
case. 

Leaving for the next section a treatment of the error involved, wc now generate data for the surviving flux at the 
91 nadir angles 90°, 89°, 0° as we did in the previous section for (80°, 40°, 0°). These data are generated using the 
complete transport equation, assuming the density given by the PREM (figure 2), and now we consider two different 
initial spectra for the sum of ATM, AGN, and GRB neutrinos (figure 6). 

Once we have the surviving neutrino flux, we have to obtain '^^}P in order to do the tomography of the Earth. The 
usual approach is to use the effective cross section to relate the surviving flux with t{9) through (7), but as we pointed 
out in the previous section, this approach is only satisfactory for very fast decreasing spectra. In this conditions, we 
expect that the reconstructed density may not be accurate in general. 

As an alternative method, we propose to use the complete integro-differential transport equation (not its solution) 
without approximation to directly obtain ^^^3^ from the surviving flux (j)i,{E, 9), which can be done noting that 

de dr ' de' ^ ^ 



As the angular variation of the number of nuclcons t is independent of the energy of the neutrino flux, if wc integrate 
over the our neutrino energy range AE = (lO'^GeV, 10''GeV) at both sides of (11) and use the transport equation (3), 
we flnd that 

de ~ D{e) ^ ^ 



where 



N{eyl 



dE 



d^,{E,e) 



Jae de 

„ danc{E',E) 



D{e)J dE( f dE' 
Jae \Je 



dE 



ME',e)-atot{E)ME,l 



Inserting -^e given by (12) into (10) yields the density reconstructed under the transport equation approach. 

The densities recovered are plotted in figure 7, together with the original one given by the PREM. Wo would like 
to point out again that wc reconstruct the earth density using the differential transport equation. We do not need to 
use the corresponding solution. As the differential transport equation is locally valid we have considered it and the 
corresponding measurable surviving flux on the earth surface. 

We clearly observe that the recovered density using the effective cross section approach is higher that the density 
that actually caused the attenuation, while with the transport equation approach the reconstructed density acceptably 
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reproduces the original one. We note that the effective cross section leads to worse results the flatter is the initial 
flux, as we saw in the previous section. This is due to the fact that NC regeneration is more important when there 
are "many" neutrinos of higher energy than a given one, that then may happen to regenerate to that lower energy. 
With a fast-decreasing spectra such as ATM alone, regeneration is not very important, and hence the cre//-approach 
also works for recovering the original density. However, in our range of energy where absorption tomography is in 
principle possible, other sources of neutrinos also are expected to contribute (AGN, GRB, etc.) and dominate in a 
total neutrino flux in that range. This justifies the use of the complete transport equation approach rather than the 
standard one. 



IV. ERROR PROPAGATION 



In this section we present the error propagation analysis considering only the uncertainty in the surviving neutrino 
flux. We have not taken into account theoretical uncertainties in the cross section such as the error due to extrapolating 
to high energies, or the uncertainty accounting for possible new physics. 

The method we suggest for performing a neutrino absorption tomography clearly depends on the surviving neutrino 
flux at each nadir angle. IceCube [11] in the southern ice seems to be the most promising neutrino detector, with 
IkmP of detection volume with photomultipliers placed regularly to catch the Cherenkov light emitted when CC 
produced muons travel in the ice. The angular resolution is expected to be of about 1*^ (or better), which is the 

angular resolution we implemented for doing the density extraction. 

As the uncertainty in a total measured neutrino flux depends among other things on the surviving flux being 
measured, we do not know what it will be for certain. Still, to test the recovery procedure, we assume three different 
uncertainties in the detected flux which are in the order of the expected ones in [11,12], and propagate each imcertainty 
to obtain the corresponding one to the recovered density. To do this, we implement the transport equation approach, 
in which is given by (12). So for an uncertainty in the flux A(p^{E,6) we find that the uncertainties in N{6) 

and D{6) are given by 



= To 



dEA(j)^{E, 



AE 



AD{0) = ^{Adi{e)Y + {Ad2{e)f 



(13) 
(14) 



where 



Adi{9)= dEatot{E)A(j},{E,T) 

J AE 

-,/ danc{E',E) 



Ad2{e) = [ dE f dE' 

J AE J E 



dE 



ame',t) 



The absolute error in 45 is then 



V D{e) J \Di9) 



(15) 



and by means of (10), we obtain the uncertainty in the number density recovered: 



An{r) = / 

^ J arcsin 



de' 



^/E? sin^ I 



dT{e') 



d6' 



(16) 



that yields the corresponding imccrtainty in the density through Ap{r) = An(r)/iV"^, which we illustrate by error 
bars in the figures 8 and 9 for the three uncertainties of 15%, 10%, and 5% in the surviving neutrino flux. We 
show the results with the error bars considering only the initial flux produced by AGN (Stecker-Salomon), GRB and 
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ATM neutrinos added together, as they are substantially the same the other initial spectra considered before (AGN 
(Protheroe) + GRB + ATM) under the transport equation approach. 

We observe that the error propagated to the recovered density is about 1.6 - 1.7 times the error in the flux detected. 



V. FINAL REMARKS 



We have presented in this work a simple method to extract the Earth's density based on the complete transport 
equation, without using any approximation as the use of the effective cross section would imply. We have compared 
the results obtained by our method with the ones obtained by the cJe//-nicthod. We can conclude that our method 
allows an accurate reconstruction for each initial flux we have considered, while the standard approach trends to 
overestimate the density accounting for its inaccuracy in dealing with the NC regeneration. 

To estimate the uncertainties in the reconstructed density, we have propagated the relative uncertainty in the 
surviving neutrino flux and found that it is increased about 60% - 70% by the procedure. 



APPENDIX A: RADON TRANSFORM AND THE EARTH'S DENSITY 



The 2-d Radon transform [13] of a continuous function n{x, y) can be defined as 

Rn{p, i)= i dx dy n{x, y) 6{p-x-£) (Al) 



Rn{p,i) (see figure 10) is thus n{x,y) integrated along the straight at a distance p from the origin and normal to 
(, = {cos 6, sin 6) which is an unitary vector defined by ^ = (cos^, sin 6), where the angle 6 is measured from the x-axis 
(Fig.3). 

If Rnip,^ is continuous, then n{x,y) can be recovered using the inverse Radon transform. The latter can be 
obtained relating the 2-d Radon transform with the 2-d Fourier tranform i^„(A^ 

Fn{Xi) = [ dxdy e'^^«n(f), (A2) 

where x = {x,y). 

Since dp 6{p - x ■ C) = 1, 

Fr,{Xi)= f dpd{p-x-i) f dx dy e'^^Q{x) (A3) 

and as the Dirac's delta 5{p — x ■ ^) allows us to write p instead of ^ • ^, we obtain 

Fn{Xi) = I dpe'^P [ dx dy 5{p-x- () n{x). (A4) 

That is, the 2-d Fourier transform of n{x) is equal to the 2-d Fourier transform of the 2-d Radon transform of n(x), 
which is known as the Fourier Slice Theorem: 

/oo 
dp e'^ P Rn{p,^). (A5) 
-OO 
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Then, by means of the inverse Fourier transform we recover n{x), 

n{x) = j^j d\K) e-^^-«X(AO (A6) 

Noting that the unitary vector ^ = (cos 0, sin cf)) , wc use the polar coordinates (A, 0) and ( A6) becomes 



1 

c^r Jo 



1-K 



"(^) - / ^ / # e-'^-«F„(AO. (A7) 







Expressing i^„(A|) according to (A5), defining p' := p — x ■ ^, and Rn{p', 1) as 

— ^ 1 /■^'^ 

we can write n{x) as 



(A8) 



-1 /"OO /"CXD 

n{x) = —J^XdXj dp e'^PRn{p,0. (A9) 

Expressing the above equation in terms of the function sign{X), performing integration by parts in p, and using 
the Fourier transform of sign{X), we obtain 

J_^ p 

which expresses n{x) in terms the p derivative of its Radon transform in a simphfied fashion. 

In the foUowing, we use the above results to relate the Earth's density with its Radon transform in a simple 
expression. 

We assume an spherically symmetric number density for the Earth which we denote by n{r) and begin by noting 
that the amount of nucleons r(z) found along a path of depth z = 2Rcos9 in the Earth is nothing but the Radon 
transform of the number density n(r) of the Earth: 

pZ poo pOO 

T{e) = / dz'n{z') = dx dy n{x) S{p - 1*- x) (All) 

Jo J —oo J —oo 

where now p = Rs'm and R is the Earth's radius. In the polar coordinates (r, <?!>) with the origin in the Earth's 
center, ^ • x = r cos(^ — then we may write 

t{9)^ d(j) dr r n{r) 5{p-rcos{4)-9)). (A12) 



Jo JO 

We can now integrate in <j) making use of the Dirac's delta property 5(/(0)) = ^^tif^^ where (f>i are the zeros of 
/((/)), and as this Dirac's delta in (A12) implies that 

p 

cos(<j)i — 9) = - then \p\ < r, 
r 

and expression (A12) becomes 
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ri9)= r^^:^kl£. (A13) 



Then as n(r) vanishes for r > R, 

^JJ2 sin^ e - sin^ 6 
We can recover n(r) by means of the expression for the inverse Radon transform (AlO): 



r2TT 



defining p' := p + ^ ■ x, we can write 



{2nfJ.^L p'-rcos{4>-e) dp' ■ ^^^^^ 



Now, the integral in (j) is defined for p' > r and p' < — r, then by the Cauchy theorem it turns out that 

( 2-n- 
d(p I p'2-r2 



2- A A, I ;y^r72 sip'>r 



p' -r cos{4> -9) — 22_ / ^ _^ 

Then, as it must be \p'\ < R so that Rn{p') does not vanish, defining 9' such that p' = -Rsin 9' , and recalling that 
Rn{p') = '''(^')' obtain that 



J arcsin 



de' drje') 

\/R^ sin^ 61' - r2 c?^'' 



(A17) 



which allows us to express the Earth's density in terms of the slope in the amount of nucleons r. 
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Figure Captions 



Figure 1: Initial muon neutrino fluxes plotted separately: AGN(Protlieroe) [7], AGN(Stecker-Salonion) [6], GRB [8] 
and ATM [5]. 

Figure 2: Earth density as predicted by the 'Preliminary Reference Earth Model' [9]. 
Figure 3: Neutrino path towards the detector. 

Figure 4: Shadowing factors for ATM and GRB fluxes for nadir angles 80", 40°, and 0°. 

Figure 5: Shadowing factors for AGN (S-Salomon) and AGN (Prothcroc) fluxes for nadir angles 80", 40°, and 0°. 

Figure 6: Total initial neutrino fluxes from AGN + GRB + ATM. AGN(P) stands for the flux by Protheroe, and 

AGN(S-S) is the neutrino flux predicted by Stecker and Salomon. 

Figure 7: Recovered densities compared with the PREM density (solid line). Boxes and stars represent the densities 
obtained with the transport equation and with the a^ff methods respectively. 

Figure 8: Densities recovered with error bars for uncertainties of 15% and 10% in the surviving flux. 
Figure 9: Density recovered with error bars for an uncertainty of 5% in the surviving flux. 

Figure 10: TZ^ region where the Radon transform is defined. The unit vector ^ and the number p define the 
integration path. 
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